Model of coherent optical spin manipulation through hot trion states in p-doped 

InAs/GaAs quantum dots 



m 
o 

(N 

On 
(N 



c3 



I 

C 

o 
o 



> 

00 

o 
o 



x 

S3 



G. Slavcheva 1 ' 2 

1 Blackett Laboratory, Imperial College London, 
Prince Consort Road, London SW7 2AZ, United Kingdom 
4 Mediterranean Institute of Fundamental Physics, Via Appia Nuova 31, 00040 Rome, Ital'Q 

A new generalised group-theoretical approach, based on master Maxwell-pseudospin equations, 
is proposed to explain recently observed enhanced circular dichroism in the excited state emission 
from p-doped quantum dot ensembles under resonant circularly polarised excitation into hot trion 
states, herein referred to as "spin-filtering effect". The theory agrees remarkably well with polarised 
time-resolved photoluminescence experiments, yielding largely unknown inter- and intra-shell spin 
relaxation time scales. This approach allows to predict optimum pulse parameters for control of spin 
dynamics, which will enable exploitation of the effect in all-optical spin-based quantum technologies. 

PACS numbers: 78.20.Bh,78.47.jd,78.47.da,78.67.Hc,71.35.Pq,71.70.Gm,72.25.Fe,72.25.Rb,42.50.Ct 



I. INTRODUCTION 



The ability to exercise strict control over the quan- 
tum dot (QD) emission properties is paramount for ap- 
plications of QDs in quantum technologies and requires 
in-depth understanding of their energy level structure. 
Three-dimensional confinement of carriers in charged 
quantum dots eliminates many bulk spin relaxation 
mechanisms but others, such as the exchange and hy- 
perfine interactions, are enhanced. Optical circularly po- 
larised excitation of charged QDs creates excited (hot) 
charged exciton (trion) states with a complex fine struc- 
ture, due to exchange interactions^. Exchange interac- 
tions play an important role in lifting the level degeneracy 
and determine the energy-level configuration and split- 
ting. The fine structure and the entailed optical selection 
rules and spin relaxation mechanisms become critical for 
successful implementation of QD-based optical spin ma- 
nipulation schemes, based on excited trion transitions^. 
On the other hand, the use of radiative cascades from 
QD excited states is a principal method for generation of 
quantum entanglement for applications of QDs as single- 
photon emitters and polarisation-entangled sources. The 
exchange interaction and the fine structure splitting of 
the energy levels involved in radiative cascades from QD 
excited states are fundamental as they determine the 
polarisation and entanglement of the emitted photons 4 . 
The fine energy level structure is revealed in the excited 
charged excitonic spectra. For these reasons interest in 
the excited charged excitons is growing. However, unlike 
the ground trion 'singlet' negative (X~) [5-8] and pos- 
itive (JT + )[9-ll] QD state that has been relatively well 
studied and understood, there have been very few reports 
on the negative hot trion, A" - *, in n-doped [12-18], and 
even fewer on the hot positive trion, X + *, in p-doped 
QDs [19-21]. 

Optical manipulation of a single electron (hole) spin, 
confined to a QD, through the resonantly driven charged 
exciton (trion) ground singlet transition, is currently con- 
sidered as one of the most promising schemes for imple- 



mentation of spin-based quantum computing, due to ex- 
tended spin lifetimes, limited by the hyperfine interaction 
[22-25]. The requirement of a resonantly driven ground 
trion transition in an inhomogeneously broadened ensem- 
ble of charged QDs, however, represents an obstacle for 
the scalability of architectures based on this quantum 
system. 

Here I show that the limitation of resonant excitation 
of the ground trion singlet state can be overcome by 
an optical excitation into the excited trion states, tak- 
ing advantage of a recently discovered effect, herein re- 
ferred to as " spin-filtering effect" . The latter consists 
of an enhanced photo-induced circular dichroism in the 
excited state emission of p-doped QDs, relative to the 
ground state, under non-resonant circularly polarised op- 
tical excitation™. The degree of polarisation is nearly 
doubled when resonant excitation into an excited dot 
level is used, due to increased spin injection efficiency. 
In addition, non-monotonic dependence of the degree of 
spin polarisation on the optical pulse power is observed, 
allowing to maximise it by optimising the pulse charac- 
teristics. This spin-filtering effect is promising for reali- 
sation of high-fidelity schemes for all-optical spin manip- 
ulation, since the increased polarisation degree contrast 
for the two pulse helicities enables highly efficient selec- 
tive excitation and readout of the spin-up and spin-down 
populations. 

In this work I theoretically demonstrate the possibility 
of using the hot trion states in p-doped QDs for effi- 
cient spin manipulation by short circularly polarised op- 
tical pulses. In Sec. II I give details of the polarised 
TRPL experiment and develop a comprehensive theory 
and a quantitative dynamical model of the optically in- 
duced circular dichroism, observed in the time-resolved 
photoluminescence (TRPL) from ensembles of p-doped 
QDs under pulsed optical excitation. I show that the en- 
hanced photo-induced circular dichroism of the excited 
state emission, compared to the ground state, stems from 
the specific trion fine structure, comprising two sets of de- 
generate states, dynamically coupled through spin flip- 
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(Color online) Experimentally detected polarised 



TRPL traces for a and a + excitation. 



flop processes, which induce asymmetry in the allowed 
optical transitions for left- and right-circularly polarised 
light. In Sec. Ill I summarize the results from the numer- 
ical simulations carried out on a realistic single QD layer 
structure and compare the calculated polarised TRPL 
traces with the experimental ones directly in the time 
domain. In particular, I discuss the dependence of the 
experimental data on detector type and the importance 
of convolving the simulation data with the detector re- 
sponse function for a proper comparison with the exper- 
imental data. The experimental data are fitted by the 
model, thereby yielding largely unknown intra- and inter- 
shell spin relaxation times, summarised in Table 1, and 
an estimate of the magnitude of the effect is provided. 
The concluding section IV summarizes the main results 
and outlines future model applications. 



II. THEORETICAL MODEL 

The QD samples employed in this work are grown by 
MBE on a semi- insulating GaAs (001) substrate 2 ^. The 
QD layers are sandwiched between two GaAs barriers 
and the dot areal density is 2 x 10 10 cm -2 with an aver- 
age uncapped height h = 4nm (N a = 5 x 10 22 m~ 3 ). PL 
experiments under resonant left- (cr~) or right- (<r + ) cir- 
cularly polarised excitation into excited dot states with 
pumping wavelength X res = 1065 nm were performed on 
QD ensembles, nominally doped with one hole. The po- 
larised PL is detected at the ground singlet X + transi- 
tion, exhibiting a peak at Xdet — 1148 nm. The TRPL 
experiments are carried out with short (r = 50 ps) opti- 
cal pulses and the PL decay in time is given in Fig. [TJ 
clearly showing marked difference in the spin polarisation 
degree for both pulse helicities. 

To understand the observed complex polarisation dy- 
namics, we consider quasi-resonant o~ or a + circularly 
polarised excitation of a p-doped QD ensemble into the p- 
shell in the presence of a resident s-shell hole, and subse- 
quent cascade relaxation to the bright trion ground (sin- 
glet) state, whose decay is detected as PL. The optical 
excitation creates X + * states, consisting of one electron- 
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FIG. 2: (Color online) Energy levels and spin configurations 
of the hot X + * trion states under a~ (left panel) and a + 
(right panel) optical excitation (upward blue wavy arrows). 
Small solid (red) and open (blue) arrowheads represent s-shell 
electron and hole total spin projection J z ; large open (blue) 
arrow - p-shell hole spin. Downward dark red wavy arrows 
denote the PL; shaded blue ellipse: exciton coupling in the 
radiative states. The excited singlet is split from the triplet 
by Ahh', the excited triplets are split by A e h- A s _ p -energy 
separation between p-shell and s-shell trion states. 



hole pair in the s-shell and a resident hole in the p-shell 
(le 1 lh 1 2h 1 ), grouped in four degenerate doublets 1 . The 
energy-level diagram and the spin configurations of the 
X + * trion is given in Fig. [3J The isotropic hole-hole 
(h — h) exchange interactions give rise to the singlet- 
triplet splitting, Ahh ~ 1— lOmeV. The isotropic part of 
the electron-hole (e — h) exchange interaction splits the 
triplet states into three equally spaced levels with spac- 
ing A e h ~ 0.1 — 1 meV. The anisotropic e — h exchange 
interaction (AEI) leads to mixing of these states. 

The system Hamiltonian of a circularly polarised pulse, 
resonantly coupled to an ensemble of 6-level resonant ab- 
sorbers with density N a , is given by: 

= (1) 
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where =p correspond to a (er + ) polarisation, fl x = p-^S 

Sl y = p^- are the time-dependent Rabi-frequencies as- 
sociated with E x and E y electric field components and 
wo and p are the resonant transition frequency and the 
optical dipole matrix element of the ground to excited 
singlet state (|1) — > |6)) transition. 

Spin decoherence mechanisms that need to be consid- 
ered include e — h AEI, which remains relevant for excited 
trion states, and the hyperfine interaction between nu- 
clei and either electrons 2 ^ or holes through dipole-dipole 
interaction 2 ^. The two degenerate 6-level systems, for a~ 
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FIG. 3: (Color online) Discrete-level model of the X + * states. 
The two degenerate six-level systems are coupled via spin 
decoherence mechanisms: 7;, - hole-spin decoherence, 7 e - 
electron-spin decoherence, 7// - spin decoherence due to spin 
flip-flop processes (AEI). Circularly polarised transitions are 
designated by O or O. Upward blue arrows: a~ x and a^ x reso- 
nant pumping, downward orange arrows: radiative or nonra- 
diative transitions; Dark red arrows: detected polarised PL. 
Curved green arrows: spin flip-flop coupling due to AEI. 



and cr + -excitation, are coupled through transverse spin 
decoherence mechanisms, as shown in Fig. [3J 

The time evolution of the 6-level quantum system in 
the presence of relaxation processes (Fig. [3]) , under time- 
dependent external perturbation (laser pulse) is governed 
by a master equation^ for the ./V 2 — 1 = 35-dimcnsional 
real state (pseudospin) vector, Sji 
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where Xj are the generators, / is the fully antisymmetric 
tensor of the structure constants of SU(6) group and 7 is 
the torque vector given in Appendix A, Eqs. (A5), (A6) 
and (A8), respectively. Tj are phenomenological non- 
uniform transverse spin decoherence times describing the 
relaxation of Sj toward their equilibrium values Sj e . The 
longitudinal spin relaxation is incorporated in the matrix 
a = diag{Tr(Y ip)) , where p is the density matrix and 
Ti,i = 1,...,6, are spin relaxation rate matrices, given 
explicitly in Appendix B, Eq. (B2) for a~ and Eq. (BIO) 
for cr + excitation. 

The excited triplet states in the two degenerate sys- 
tems, where the photo-generated and the resident hole 
spins are parallel (Fig. |2J), are coupled via spin flip-flop 
processes, whereby the electron and the hole flip simulta- 
neously their spin (level |3) and |4) in Fig. [3]). This e — h 
AEI-mediated process plays a significant role in n-doped 
and p-doped QD spin dynamics 2 ^. Due to dynamical cou- 
pling of the two degenerate systems through spin flip-flop 
processes and simultaneous electron and hole spin rever- 
sal, the total spin projection of these states changes and 
as a result additional transitions become allowed due to 
spin selection rules. For instance, the total spin projec- 
tion of state 1 3) changes from | — 7/2) — > | + 1/2), and of 
state 1 4) from | — 5/2) — > \ — 1/2), thus making possible the 
transitions |6) |4), |6) |3), |4) |2) and |3) -> |2), 
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FIG. 4: (Color online) Energy-level diagram for a + excitation, 
modified by the spin flip-flop processes in the triplet states 
(level 1 3) and |4), right panel). Downward (orange) arrows 
denote radiative and non-radiative transitions before spin-flip 
flop. The green double-headed arrows represent the spin flip- 
flop of the electron and the hole, changing the total spin of 
state |3) from | - 7/2) -> | + 1/2), and of state [4) from | - 
5/2) — > | — 1/2). This results in additional allowed transitions, 
denoted by red arrows. 
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FIG. 5: (Color online) Refractive index profile of the simula- 
tion domain, containing InAs QD layer (h = 4nm), embedded 
between two GaAs barrier regions, 50 nm each. The circularly 
polarised pulse, modelled by two x- and y- linearly polarised 
waves, phase shifted by 7r/2, and its E-field components, are 
schematically shown. 



denoted by red arrows in Fig. SI Therefore more transi- 
tions at any one time become allowed in one polarisation 
(here without loss of generality chosen as a + ) than in the 
other (er~). 

The simulations domain consists of an InAs QD layer 
with nominal thickness given by the height of the typical 
QD [h = 4nm), embedded between two GaAs barrier 
regions each with thickness 50 nm (see Fig. [5]). 

Maxwell's curl equations for a circularly polarised laser 
pulse propagating along the QD growth direction, chosen 
as z- axis read: 
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where P x and P y are the macroscopic polarisation vector 
components. Eqs.fl2Jl and ^ are coupled through the 
macroscopic polarisation induced in the medium by the 
electromagnetic wave (see Appendix C): 
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where =F corresponds to cr" (<r + ) excitation. 

The circularly polarized Gaussian pulse at the left 
boundary of the simulation domain, z = 0, is modelled 
by two x- and y-linearly polarised waves with resonant 
carrier frequency, ojq, phase shifted by ir/2: 
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(5) 



where E$ is the initial field amplitude. The system 
of equations ([5][5]) is discretised in space and time ( 
Az = 1 A, At = 3.33 x 10~ 4 fs) and solved numerically 
directly in the time domain using the Finite-Difference 
Time-Domain (FDTD) technique. 



III. NUMERICAL SIMULATIONS AND 
COMPARISON WITH EXPERIMENT 

A. Simulated dynamics 

An optical pulse with pulse duration r = 50 ps (se- 
lected to match the experiment) is injected at z = 0. 
The pulse center frequency ujq, is tuned in resonance 
with the energy splitting between the ground heavy-hole 
level |1) and the excited singlet trion state |6): AEi^ 6 = 
Ex+ + A s _ p + Ahh, where the ground trion singlet en- 
ergy E x + = 1.0815 eV is determined from the resonant 
PL spectra at Xdet- A s _ p ~ 73 meV is inferred from 
the PL spectra^ 6 and A^h = 12 meV, A e h ~ 0.5 meV are 
taken in agreement within. 

The pulse area below the pulse envelope is chosen to 
be 7r, corresponding to E — 2.69 x 10 5 Vm" 1 (Eq.©), 
assuming dipole matrix element p — 9.83 x 10" 29 Cm. 
The pulse power, corresponding to a 7r— pulse is equiv- 
alent to 754 mW (~ half of the peak pulse power deliv- 
ered by the Thsapphire laser). The dipole moment of 
the X + ground singlet transition, p ~ (7 ± 2) x 10 -29 , 
[30-32] is taken as initial approximation. Radiative de- 
cay rate T21 ~ 1.27 ns" 1 of the X + ground singlet 
transition is measured ir* 2 ^. Spontaneous emission rates 
T4.1 ~ 1.35 ns —1 , rsi ~ ~ 1.2 ns" 1 are estimated, 
using the Fermi golden rule^. 

Non-radiative spin relaxation times ~ 1 — 7ps _1 are 
calculated in^i, although conflicting values for the p —> s 
intcrshell relaxation times can be found in the literature: 
from fast, 1.5 — 9ps _1 to very slow 750 ps -1 — 7.7 ns" 1 . 








lb) 



500 10GO 1500 2000 250D 3OD0 3500 









m 



5M 1O00 1500 2000 2500 MM 3500 
Time [ps] 

FIG. 6: (Color online) Simulated short-time (a snapshot 
at t = 120 ps - left panel) and long-time (a snapshot at 
t — 3.5 ns - right panel) spin population dynamics of all six 
levels, following cr" (a, b) and a + (c, d) circularly polarised 
Gaussian 7r-pulse with pulse duration T p = 50 ps and E-field 
components E x and E y . The spin population of level |2), P22 
(black curve), is proportional to the detected PL emission. 



Our simulations show that slow intershell relaxation 
times do not reproduce the observed TRPL time evo- 
lution, therefore we assume T52 ~ 5 ps" 1 . The intrashell 
(s* — > s) rate Yq2 ~ 5ps" : is taken in agreement with2&. 
The timescale of the intershell (s* —> p) hole spin re- 
laxation is largely unknown. Our simulations show that 
these processes occur on a tens of ps timescale, but faster 
than 25 ps" 1 rates lead to unphysical negative spin popu- 
lation. We assume r 6 5 ~ 50 ps" 1 , correctly reproducing 
the experimental TRPL trace. 

The spin decoherence rates for the transverse relax- 
ation processes (Fig. \S§ are taken 7 e = 500 ps" 1 for elec- 
tron spin decoherence through hyperfine interaction with 
the lattice ion spins 25 , jh — 14ns" 1 for the hole spin 
decoherence due to dipole-dipole interaction 28 ; the spin- 
flip rate, Jff = iaei ~ 125 ps" 1 (AEI ~ a few tens of 
peVffi). 

For a~ excitation, the dot density, N a , and the dipole 
moment matrix element, p, have been identified as key 
parameters significantly affecting the time trace. The 
non-radiative decay times, tq2 and T52, are critical for 
obtaining the correct amplitude. For a + excitation vari- 
ation of the radiative decay times, T51 and T41, signifi- 
cantly reshapes the TRPL trace. 

The set of parameters producing best fit to the exper- 
imental data is given in Table I. 

The simulated spin population dynamics at a point 
in the QD layer using the best fit parameters (Table I), 
is plotted in Fig. [BJ showing significantly different long- 
time a~ and a + TRPL decays (/?22 b, d). For com- 
parison with the TRPL data, the computed TRPL time 
traces (P22 in Fig. [5]) are plotted along with the exper- 
imentally detected photoluminescence by two detectors: 
InP/InGaAsP photomultiplier(PMT) and a microchan- 
nel plate detector (MCP) with extended SI photocathode 
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TABLE I: Characteristic parameters used in the simulations of the TRPL following &T P — 50 ps pulse excitation with a~ and 
a + polarised light. The decay timescales between a pair of levels \i) and \ j) are denoted by ry, and taei is the decay timescale 
associated with the anisotropic e-h exchange interaction; meaning of all parameters explained in the text. 
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FIG. 7: (Color online) Circular polarisation dynamics of the 
PL, following resonant excitation into p-shell states, experi- 
mentally detected from the QD X + ground singlet trion state 
by PMT and micro-channel plate (MCP) detectors, and the- 
oretically computed polarised TRPL trace for a~ and a + ex- 
citation (p22 in Fig. EJ) 



in Fig. It becomes apparent that the detected TRPL 
trace shape strongly depends on detector characteristics: 
the better the detector is, the closer it is to the simulated 
curve. 

In order to properly compare our theory with the ex- 
perimental data, we convolve our simulations with the 
photomultiplier detector response function (Figs. [T21 EH 
Appendix D). 



for the first time that the entire duration of the TRPL 
signal is crucial for correct reproduction of the PL time 
decay shape. 

The influence of the dot density on the TRPL time 
evolution trace is given in Fig. [5] for four simulations 
with the rest of parameters kept fixed. We start with an 
estimate of the resonantly excited dot density on the or- 
der of (0.2 — 0.3)^ Fig. j9ja) and increase it stepwise to 
the density per unit volume of 5 x 10 22 m~ 3 , calculated 
on the basis of the areal dot density of 2 x 10 10 cm~ 2 
and average height, h = 4nm, of the QD layer. The 
departure of the theoretical curve from the experimental 
one, after the roll-over in Fig. [S^a) is quite substantial, 
and in order to save computational time, we interrupt 
the simulations as soon as the general trend is apparent 
(the kink in the theoretical curve is due to earlier inter- 
ruption of the simulations). As the dot density increases 
by nearly three orders of magnitude (Fig. G^b)), the dis- 
crepancy becomes smaller and the agreement between 
theory and experiment is significantly improved in Fig. 
injc) and (d). The best agreement between theory and 
experiment is achieved for the calculated volume density 
N a = 5 x 10 22 m -3 . Further increase of the dot density 
to N a = 7 x 10 22 m~ 3 and beyond, does not lead to any 
improvement of the fit. 



C. Comparison with experiment 



B. Simulations length in time and dot density 
influence on the TRPL trace shape 

The calculated TRPL shape, and therefore the agree- 
ment between the theory and experiment across the en- 
tire PL decay time, is very sensitive to the length in time 
of the simulations. In Fig. [SJ the spin population, P22, 
time evolution for three different simulation lengths in 
time and the corresponding comparison between theory 
(after convolution with the detector response function) 
and experiment for a cr~ pulse excitation are shown. To 
the best of our knowledge there have been no attempts to 
look at the TRPL decay beyond 2.5 nsi 2 ^^!. We show 



The convolved theoretically calculated polarised TRPL 
traces are plotted along with the experimental data for 
both pulse helicities in Fig. [TU] for the parameters listed 
in Table I. The theoretical results agree remarkably well 
with the experiment. The rising slope and the decay im- 
mediately after the peak are almost perfectly reproduced. 
The small discrepancy at longer times can be attributed 
to several factors: there is no way, in the experiment, to 
determine precisely which hot trion state is resonantly 
excited, so the assumed X res is only approximate; A s _ p 
is inferred from the resonant PL spectra, which is only 
a rough estimate, as the excited state can be viewed as 
an envelope of a number of spectral lines; the inhomoge- 
neous broadening is not taken into account. The lower 
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FIG. 8: (Color online) Time evolution of the spin population of level |2), P22, following a~ pulse excitation, for three different 
simulation lengths in time (a) 1.36 ns; (b) 1.95 ns; (c) 3.34 ns (upper row) and the corresponding comparison between theory 
(after convolution of the raw simulation data, plotted on the upper row, with the detector response function from Fig. 4) (blue 
curve) and experimental TRPL trace (red curve) (lower row). The simulation parameters are as follows:r2i = 1.27ns, T41 = 
1.35 ns, T51 = 1.2 ns, r 52 = 5ps, r 6 i = 1.2ns,r 6 2 = 5ps, t 6 s = 50 ps, taei = 125ps,A res = 1127nm,_E Js:+ = 1.08 eV, pi->6 = 
8 x 10 _29 ,iV a = 5 x 10 22 m _3 ,£o = 3.3 x 10 5 Vm" 1 




FIG. 9: (Color online) Comparison of the simulation data after convolution (blue line) with experimentally detected TRPL 
time trace (red line) for a~ pulse excitation with pulse duration T p = 50 ps and dot density (a) N a = 5 x 10 19 m~ 3 ; (b) 
N a = 1.5 x 10 22 m~ 3 ; (c) N a = 3 x 10 22 m~ 3 ; (d) iV a = 5 x 10 22 nT 3 , with the rest of simulation parameters as in Fig. [8] 



peak height for a + excitation can be attributed to the 
estimated value for the spin flip-flop coupling rate and 
requires further investigation. This is, however, beyond 
the scope of this work. The theory correctly predicts cir- 
cular dichroism in the polarised TRPL emission for both 
pulse helicities, shown Fig. QTJ 



IV. CONCLUSIONS 

The proposed new theory and model explain quanti- 
tatively the origin of the spin-filtering effect. The ob- 
served enhanced time-resolved circular dichroism in the 
excited trion state emission is shown to emerge from the 
dynamical spin flip-flop coupling between the two de- 
generate sets of excited charged excitonic states. Re- 
markably good agreement between theory and experi- 



ment is achieved, thereby yielding largely unknown in- 
tra and inter-shell spin-relaxation timescales. The ad- 
vantageous to this effect increased optical spin injection 
efficiency under quasi- resonant excitation into QD hot 
trion states, compared to the ground trion singlet, opens 
up new avenues for high-fidelity optical spin manipula- 
tion. The model represents a tool for investigation and 
identification of an optimum pulse power for achieving 
maximum degree of spin polarisation, thus maximizing 
the efficiency of the optical spin orientation. This en- 
ables high-fidelity quantum computing gate operations. 
The proposed formalism can be extended to describe the 
spin dynamics in multiply-charged QDs for which further 
enhancement of the circular dichroism is reported^. 



7 







— cr' - theory 


s 




n" - expefimcntt 


7 

j 






5 

! J 




(a) 







90 95 

lime [ns] 




90 95 

lime [nsl 



FIG. 10: (Color online) Comparison between simulated TRPL trace (after convolution with detector response function) and 
experimentally detected TRPL trace for (a) a~ excitation and (b) a + excitation. 
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FIG. 11: (Color online) Theoretically computed polarised 
TRPL traces for a~ and <r + excitation; simulation param- 
eters from Table I, exhibiting circular dichroism. 
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Appendix A: Generators of SU(6) Lie algebra and 
pseudospin equations 

For an N- level system there are N 2 — 1 generators, Xj, 
satisfying the commutation relations: 



[Xj, Afe] = 2i fjkiXi 



(Al) 



where fjki is the fully antisymmetric tensor of the struc- 
ture constants of SU(N) Lie algebra. / guarantees con- 
stant length of the real state vector, |S| = cons&. 
For the 6-level system considered, there are 35 lambda- 
generators. A possible choice satisfying the othogonality 



relations: 



TriXi.Xu) =26 



Ojk 



(A2) 



can be constructed from the transition-projection opera- 
tors: 

P mn = |ro) (n| (A3) 
where \m) and |n) are eigenstates of the system Hamil- 
tonian, using the following definitions: 



lAjk Pjk ~t~ P kj •> 

Vjk = -i (pjk - Pkj 

2 



l<j<k<N 



Wl 



1/2 



Pll 



. + P n - lPi+i,i+i) ,1<1<N 
(A4) 

5. The vector A, defined as 



with N = 6 and I = 1,2 
an ordered array of the elements of u, v and w: 



X = (ui 



.,Vi2, ...,w 1 ,...,w N -i) 



(A5) 



has operator components, Xj, that satisfy the commu- 
tation relations Eq. (|Aip and are traceless. The latter 
represent the generators of the SU(6) group. 
The fully anstisymmetric tensor, given by: 



fijk = \^ [Tr [Xi.X k .Xj ) - Tr ( X l .X J .Xk 



(A6) 



has 750 non- vanishing components, and the torque vector 
is given by: 



7i (*) 



\Tr(H(t)h 



(A7) 



where the time-dependent Hamiltonian, H, is given by 
Eq. (1). For cr~ excitation Eq. (|A7Jl gives: 



7 = (o, o, o, o,-o x ,o, o, o, o, o, o, o, o, o, o, o, o, o, o,-n V) o, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

(A8) 

a a i -2A eh -A feh +A ii p+wo A cft -A hh +A S p+^o 5A e(l -A hfe +A sp +w 4A hh +A Bp +u \ 

-t\ hh - A sp +W , -j. , -j. , -7- , ^_ ) 

I 



Substituting the torque vector and the non- vanishing ten- 



sor components into the equation of motion Eq.(2) for the 
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real pseudospin vector, S: 



1,2, ...,35 



(A9) 



3 = 1 fe=l 



r 3 = 



(0 o 
1 

I II 

) 
) 
3 



II 

I) I) 

734 

74: 

I) I) 

I) I) 



0\ 

(I 





(I 

0/ 



/o o 









\0 I) 







-(r«+7i3+7«) 


II 



a set of 35 equations for the real state vector compo- 
nents is obtained, which can be written in a matrix form: 



S = MS 



(A10) 



where the dot denotes time derivative and the matrix 
M is a 35 x 35 matrix, depending on the energy level 
splittings and the Rabi frequencies, Q x and Ct y . Due to 
the problem symmetry, the pseudospin equations remain 
valid for a + excitation, provided that Q y is replaced by 

— ily. 

The dissipation in the system, through spin relaxation 
and decohcrence processes, is taken into account within 
the master pseudospin formalism^ and is explained in 
the following section. 



/ 





II 


1) \ 


/ 





\ 







































II 







745 








f 6 = 





II 







o -(r 6 . 


+ r M 


+ 754) 












\0 








r 65 / 


V o o o o o - (r 6 i 


+ r 62 + r 6 . 


)/ 



(B2) 



Then we can define a = diag(Tr(Tip)), thus isolating 
the longitudinal spin relaxation from the transverse one 
and write the equation of motion for the density matrix: 



dp 
dt 



p,H 



■<r -T t p 



(B3) 



Appendix B: Longitudinal and transverse spin 
relaxation matrices 



Consider a~ optical excitation of the |1) — > |6) tran- 
sition (Fig. 2, left panel) and subsequent relaxation via 
all dipole-allowed transitions, with corresponding decay 
rates, to the ground state with spin-up heavy hole. The 
population of each level satisfies the following set of mas- 
ter equations for the diagonal density matrix compo- 
nents: 



P, H | + r 6 i/5 66 + r 5 i/5 5 5 + T41/944 + r 2 l/5 2 2 



dpn 


i 


dt 


~ h 


dp22 


i 


dt 


~ n 


df>33 


i 


dt 


~ n 


dpu 


i 


dt 


~ Ji 


dp55 


i 


dt 


~ h 


dpm 


i 


dt 


~ h 



p,h\ + r 62 p66 + r 52 P55 - r 2 ip 2 2 

22 



= T\P' H \ + 743 P44 - 734P33 



33 



(Bl) 

(745 + 743 + r 4 l) P44 + 734P33 



Pi H I rg + r 65P66 + 745P44 ~ (754 + T 5 i + F 52 ) P55 



P,H\ +754P55 



(if) 



(^61 + r 6 2 + r 65 ) p 66 



where I\j and 7,j are the relaxation rate between lev- 
els i and j and I have introduced, for the sake of gener- 
ality, different rates of spin population transfer, due to 
e-h anisotropic exchange interaction, between the excited 
triplet levels (denoted by ^aei in Fig. 3). For each of 
the above equations we can introduce a 6 x 6 matrix, 
Ti, i = 1,2, ..,6, describing the spin relaxation rates as 
follows: 



r, - 



1 o 

) r 21 o o o 

) o o o o 

) o o r, u o 

) o o o r 5 i 

1 o 



\ 






W 



r 2 = 



o 
ii 
ii 



\o 



I) 

-r,i o o o 



I) 

o r 53 





\ 






w 



The transverse spin relaxation (decoherence) matrix, r t , 
is obtained by writing down master equations for the off- 
diagonal density matrix components, taking into account 
the spin decoherence mechanisms relevant for each level 
(see Fig. 3, double-headed blue arrows coupling the two 
6- level systems): e.g. for the ground level - hole-spin 
decoherence rate 7^, due to dipole-dipole hyperfme in- 
teraction with the lattice ion spins, for the levels with 
unpaired electron spin - 7 e , electron spin decoherence 
via hyperfme interaction, and for the excited triplet lev- 
els - 7//, spin flip-flop processes, whereby the electron 
and the hole simultaneously flip their spin. As a result, 
the transverse spin decoherence matrix reads: 



r ( = 



/ 7 e 7/j 7/, 7/j 7ft \ 

7e 7 e 7 e 7 e J e 

Iff Iff Iff Iff Iff 

Iff Iff Iff Iff Iff 

7e 7e 7e le 7 e 

V le 7e 7e le le / 



(B4) 



The density matrix, p can be expressed in terms of the 
lambda-generators of SU(N) Lie algebra, where N = 6 is 
the number of discrete levels of the quantum system, as 
follows: 



N -1 



Sj (t)A i 



(B5) 



where Sj are the real state vector components and I is 
the unit matrix. The elements of the diagonal a matrix in 
terms of the real state vector components can be obtained 
by taking the trace of the product of p from Eq. (|B5p with 
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the Ti matrices from Eq. (|B2 



of the spin decoherence rates by: 



o-ii = ^[(lO + 30S 3 i - 1()V3S 32 - 5VTjS 33 - 3VIOS34 - 27l5S ;i5 ) r 21 

+i5v / Ss 3 3r 4 i - 3VTo (r 4 i - 4r 5 i) s m - 2VT5 (r 4 i + r 51 - sr 61 ) s 3l 
+io(r 4 i + r 5 i + r 61 )] 

622 = ^,[(-10 - SOSsi + loVaSsa + 5V6S33 + 3VT0S 3 4 + 2\/l5S 35 ) r 2 i 

+2 (5 + 6\/ios 34 - VisSss) r 52 +10 (1 + t/TES3s) r 62 ] 

*33 = 7^[f-10 - 20V35 32 + 5V65 33 + SVWSu + 2^15S ;i 5) 734 
60 V J 

+ 110 + 15\/6533 - 3VT0S 34 - 2VT5S 35 J 743] 

*44 = ^r[f 10 + 2073S 32 - 5V6Sf S3 - 3VIOS34 ~ 27I5S35) 734 
60 V I 

+2 (5 + 6VTOS3.1 - VT5S35J 754 

- (lO + I5V6S33 - 3VT05 34 - 2VT55 35 ) (7-13 + 745 + r 44 )] 



^55 



= (lO + 15V^5 33 - 3^10534 - 2VT5S35) 745 

-2 (5 + sVwsu - vTs&s) (754 + r 8 i + r 52 ) 
+10 (1 + v^5S 35 ) r 65 ] 



5-66 = - g (1 + VEs 35 ) (pei + r 62 + r 65 ) 



(B6) 



Substituting a into the second line of Eq.(2), we can 
find expressions for the relaxation times of the spin pop- 
ulation terms towards their equilibrium values, which for 
<j~ excitation read: 



T\ — Tg — T 15 — T 16 

T2 = 13 = T17 = r lg 



Ty-i — Toa — T;in — 



T A = T, 



■T 2t 



2 



2 

6 — J 7 — 1 & — J 21— J-22 — - 



(B8) 



T10 



T 25 



1 

W 



: T14 — T26 — T27 — T21 



To: 



In the case of a a + excitation, the energy-level diagram 
is more complex, as additional relaxation channels are 
dynamically switched on (Fig. 4). The population of 
each level satisfies a different set of master equations for 
the diagonal density matrix components: 



dpn 
dt 

dp22 

dt 
dp-i-3 
dt 

dt 

dt 
dpm 
dt 



h v ' H \ 11 + r61l ° 6B + rsipss + r 4 ip 44 + r 2 ip 2 2 



- rfi2P66 + T52P55 + ^42pM + f 3 2P33 ~ T21P22 



= j [p, fl] ^ + r 63 p 66 + r 43 p 44 - (734 + r 32 ) p 33 (B9) 

= ^ [P> H\ 44 + r 64PB6 + 754P55 - (745 + 743 + T 4 i + T 42 ) p 44 + 734P 33 



p, h\ + r 65 p Be - (754 + r 5 i + r 52 ) p 55 + 74sp 44 



= P 



(Psb + T64 + 1^63 + ^62 + Tei) P66 



We can introduce, similar to the a~ case, 6x6 matrices, 
Fj, i = 1,2,.., 6, describing the spin relaxation rates as 
follows: 



T31 
I32 
T33 
T34 
T 35 



2 

3 

734 



24 



4 734 + 12 743 + 9 745 + i2r 4 i 

8 

745 + 4 (754 + r 5 i + r 52 ) 
2 



r 6 i + r 62 + r, 



(>5 



(B7) 



Ti = 











(i 





\ 
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-r 2 i 
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fa = 





r 32 




















r« 













r 12 
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r 5 i (i 















r 62 
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r 6 
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Vo 


n 


11 
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(B10) 
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fj - 


o 


734 



















7 13 
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- (74s + 743 + Tn +T4 2 ) 
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r 63 y 




V o 











r 6 n 
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/ n 





\ 
















n o 





















o 










745 








fu - 


n o 










-( 


r 5 i + r 52 


+ 754) 
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\ n 








r K / 


\ — (r,ii 


+ r ea + r B3 + r 


84 + r e ,,)/ 



The transverse spin decoherence times, Tj,j = 
1, 2, 30 in the first line of Eq. (2), are given in terms 



The longitudinal spin relaxation can be incorporated in 
a diagonal matrix, <r, whose components are calculated 
as the Tr(Tip) and subsequently expressed in terms of 
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the real state vector components, Sj,j = 31, 32, 35, 
responsible for the spin population difference between the 
pairs of levels. The longitudinal spin relaxation times in 
Eq.(2) are given by: 



T 3 i 
T33 
T 35 



2 



3r 



.32 



2 7 ; 



.",4 



24 



4 734 + 12 743 + 9 745 + 12 (r 41 + r 42 ) 

8 

745 + 47 54 + 4(r 51 + r 52 ) 
2 



(Bll) 



r fi 



r 



62 



r 6 4 



The transverse decoherence times are given by 
Eq.flBi. 



Appendix C: Relationship between the macroscopic 
polarisation and real state vector components 

The macroscopic medium polarisation is given by the 
expectation value of the dipole moment operator, e(Q), 
where Q is the local displacement operator. 



P = -N a e (Q) = -eN a Tr(p.Q 



(CI) 



where N a is the density of the polarisable atoms (QDs) 
in the medium. Since the electric field vector of a cir- 
cularly polarised pulse is in a plane perpendicular to the 
propagation (growth) z-direction: 



E = E T 



Ey e y 



(C2) 



the dipole moment induced in the medium will be along 
E x and E y electric field vector components, therefore the 
local displacement operator can be decomposed in two 
components, along x- and y- axis, selecting appropriate 
lambda-generators in such a manner that they produced 
a coherent, circularly polarised excitation term in the 
dipole interaction Hamiltonian: 



Q = Qo ' 



/ 1 \ 




















> 























































e x + 




















e y 



























\ 1 / 




V i 














/ 





3o 



(A 5 e x + X 2 oe y 



(C3) 



where e x and e y are the unit vectors along x— and y— 
axes, qo is the typical distance between the charges in a 



dipole and the dipole moment is given by p = eqo. With 
this choice of lambda-generators and substituting E and 
Q from Ea. l|C2p and (|C3|) . the dipole-coupling interaction 
Hamiltonian acquires its form from Eq.(l): 



flint (t) = -eE.Q = -eq [E X X 5 



Ey A20 



= h 



\-Un x ±in y ) 



/ 



. (C4) 

where the upper (lower) sign denotes a (er + ) excitation 
and we have defined the complex Rabi frequency compo- 
nents fl x = p-fs fly = p-^. 

Taking the trace of the product of the density ma- 
trix operator in terms of the real state vector, given by 
Eq.(B5) and the local displacement operator Eq. (C3), 
one can obtain Eq.(4). 



Appendix D: Convolution of the simulation data 
with detector response function 

The detector response function of the InP/InGaAsP 
PMT detector, as given by its specifications, shown in 
Fig. [T3] (a), is converted in the time domain, Fig. [T5] 
(b). Since the detector sampling rate Ati = 0.05 ns 
is much coarser than the simulations sampling rate of 
At2 = 3.33 x 10~ 4 fs, we add n = [Aii/A^] points within 
each bin and thus discretise it, Fig. [T2Tch 

In order to calculate the convolution of the simulation 
data, (?, with the detector response function, h: 



g*h= J g(r)h(t-T)dT (Dl) 

— OO 

we take advantage of the Convolution Theorem: 

F(g*h)^G(f)H(f) (D2) 

where G and H are the Fourier-transforms of the con- 
volved functions, respectively g, and h, and calculate 
inverse Fourier transform of the product of Fourier- 
transforms above. The simulation data is padded with 
zeros to match the length (in number of points) of the 
detector response function and is convolved with the lat- 
ter. The functions before and after the convolution are 
shown in Fig. [T3J 
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Number of lime bin time [ns| tirrw |ns] 



FIG. 12: (Color online) (a) InP/InGaAsP PMT detector response function, as given by the specifications; (b) Converted in 
the time domain with sampling rate Ati = 0.05 ns; (c) Discretised with step size of the simulation sampling rate, Ai2 = 
3.33xl0" 4 fs. 




iime[ns] time [ns] time(r«] 



FIG. 13: (Color online) (a) Simulation data for P22 versus time; (b) Simulation data padded with zeros to a length in time 
matching the one of the detector response function; (c) Convolution of the simulation data in (b) with the discretised detector 
response function from Fig. I12f c'l. Note that as a result of the convolution, the signal is translated along the time axis. 
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